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Motivation: Learn More About 
lonospherically Induced 
Frequency Variations 


On-Air Spectrograms of 5 MHz WWV 
Showed Smooth Daytime and 
Turbulent Nighttime Characteristics 
with a Prominent Negative Swing at 
Dusk and Positive Swing at Dawn. 
The Dawn Swing Additionally 
Showed Multiple Higher Order 
Swings, that Sometimes Manifesting 
Abruptly Mid-transition 


Mode Splitting Observed on WWYV During 
Dawn and Dusk Transitions 


Positive Frequency Excursions During Sunrise 


Ray Trace Simulation Predicting 
Multiple Propagation Modes from 
WWV to WASFRF After Sunup 


Postulated mode splitting mechanisms: 


1. For multiple modes returning from a common 
layer, total path length increases according to 
number of hops. The longer modes close faster 
for a given layer descent speed, resulting in more 
Doppler shift according to number of hops. 


2. Multiple hop modes that manifest abruptly do 
so when the required high angle of propagation is 
supported. 


3. Bill Engelke AB4EJ suggested different Doppler 
tracks for simultaneous refractions from different 
layers. 


4. Bill Liles WQ6Z suggested changing wave speed 
can be responsible for producing frequency shifts. 


2D Ray Trace from WWV to WASFRF, 4-50° in 1° steps, January 2020 at 1430 UTCon 5 MHz 
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Simulation Provided by Carl Luetzelschwab K9LA 


Measurement Tools and Techniques to Unravel 
Multiple Processes by Specific Testing 


Method to obtain high resolution frequency data: offset receiver low in frequency in 
USB mode to down convert the RF carrier to an audio tone. This preserves small 
frequency variations on that can be analyzed on a sound card spectrum analyzer or an 
audio frequency demodulator. A GPSDO provides atomic clock accuracy and stability. 
Compare simultaneous WWV spectrograms across multiple frequencies. 

Compare oscillographic time domain presentations of a GPSDO stabilized tone to the 
down-converted WWYV carrier for study of short term amplitude and frequency 
variations. 

Measure Times of Flight of multiple modes using WWV per-second timing ticks and 
GPSDO synchronization to deduce propagation modes. Correlate data with ray trace 
simulations using actual solar-terrestrial data. 

Separate simultaneous WWV and WWVH spectrogram data by deinterlacing the 500 
and 600 Hz side tone encoding. 


Comparison Between Spectrogram and Single Frequency 
Formats During HamSCl Festival of Frequency Measurements 
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Single Frequency Detector using FLDIGI 
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Path: WWYV, Ft. Collins CO to WASFRF, near San Antonio, TX 0000z — 2359z October 1, 2019 6 


Simultaneous WWV 
Spectrograms at 
2.5,5, and 10 MHz 


During Dawn Transition 


2.5 and 5 MHz frequency steady during #7 
the day and turbulent at night. 10 MHz _ | : 
frequency turbulent both day and night. (Sede 
Only weak correlation over frequency pee CO-1306z 
during turbulence. The turbulence sow V%.12322 
amplitude does not scale with carrier | : 
frequency. 


2.5 and 5 MHz manifest multiple high 
order modes during sunrise but 10 MHz 
does not. 


Some modes manifest or extinguish 
abruptly mid-transition. 


Ft. Collins, CO (40.68, -105.04) to San Antonio, TX (29.57, -98.89) on October 11, 2019 


All f 
requencies show a simultaneous GPSDO on 2.5 MHz receiver only. 


mode showing little frequency shift. 


Time Domain Video Displays for Analysis of Short Term Variations 
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Rapid in-phase, null, and out-of-phase sequence occurring in less than 100mS 


Two 1 kHz tones with atomic accuracy and stability: A GPSDO stabilized 
function generator and on-air 5 MHz WWV downconverted to 1 kHz and 
filtered to 10 Hz bandwidth. The difference is that WWYV arrived by skywave. 


This type of comparison in the time domain succinctly demonstrates AM and 


FM modulations imposed by the ionosphere during skywave propagation. Spectrogram Showing 


Observed behavior includes steady amplitude with little frequency variation, symmetrical Frequency 


rapid amplitude flutter with frequency shifts on the order of Hertz, and even 
larger amplitude fluctuations that include complete carrier nulls with phase 
reversals. 


Spread Suggestive of 
Modulation Sidebands 


Rapid and deep amplitude fluctuations can be considered a modulation 
process that can superimpose modulation sidebands on spectrograms. P 


Passing Downconverted Carrier Through a Hard Limiter 
Sharpened Spectrogram by Stripping Amplitude Information 
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Time of Flight Measurements Using WWV 
Per-Second Timing Ticks 


Published Timing Tick Information: 
WWV -—5cycles of 1000 Hz 
WWYVH — 6 cycles of 1200 Hz 

Alternating 500/600 tone frequencies 


A sync reference can be obtained from 
the 1 pps output from a GPSDO. 
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On-the-air reception of WWV and WWVH Timing Ticks 
At WASFRF. Time difference of arrival is ~17mS. 
WWYV at 1350 km, WWVH at 6050 km great circle dist. 
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Idealized Geometry for 1, 2, and 3 Hop Paths Used to 
Approximate Expected Pulse Timing 


Geometric Times of Flight for 1, 2, and 3 Hop Paths from 


— he=virtual height WWYV to WASFRF vs. Reflection Layer Height 


P, = 2\/(X/2)2 + h2 
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Reflection Layer Height - km 


For ground distance (X) of 1350km and reflection 
heights (h) from 150-300km, the differences in arrival 
times between the 1 and the 2 & 3 hop modes range 
from 0.3 to 4 milliseconds. This means the 2-hop and 3- 
hop pulses will come down within the 5 mS Primary 
pulse length, resulting in superposition. 


Simplifications: 

1. “Flat Earth” 2D geometry 

2. Perfect reflections at constant 
virtual height 

3. Constant wave speed =c 
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The Time Delays Between Multiple Modes are Less Than the Pulse Width. Therefore the 
Primary and Delayed Pulses Overlap, Lengthening the Pulse by Superposition. 


Pulse 1 


Simulated |i a oe or ee 


1 mS Delay -> 6 —_ 2 mS Delay -> 7 a 3 mS Delay -> 8 cycles 
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Pulse Extension 


Measured © 
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Timing jas was extracted _ measuring time eae from the GPSDO sync to the miesdine edge for the 
primary pulse and by the pulse extension on the trailing edge for the delayed pulses. Timing data was 
acquired over the 4 hour dawn transition on January 29, 2020. 
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UTC on January 29, 2020 


- Negative slopes are consistent with Faget ere EPSwolaync coe name eae) 
decreasing path lengths from descending Corrected for 4.75mS Receiver DSP Delay 

ionization layer. 

- Abrupt appearance of 3 Hop mode is 
consistent with spectrogram data and 


high angle propagation theory Pas) eta 1 Hop 
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Actual arrival times of second and third delayed pulses can be estimated by adding the 
delay times in the top graph to the Primary arrival times in the bottom graph. 
Effective virtual height descent rate available from slopes. 13 


Measured TOF’s were Consistent with Theoretical TOF’s for 
a 255km Refraction Layer and Support Premise that 
Delayed Pulses Are from Multiple Hops 


Geometric Time of Flight for 1, 2, 3, and 4 Hop Paths from Measured and Interpolated TOF's from GPSDO Sync to 
WWYV to WASFRF vs. Reflection Layer Height Primary and Multiple Modes 
g Matched , 
4 Hop Scales WASFRF Sunrise 1325z | | WWYV Sunrise 1414z 


7 
3 Hop 
6.5 7 


EEE 


-mS 


| 
a 
fe) 
SG 
wo on 
in 
Geometric Time of Flight - mS 
Measured Time Of Flight 


285 275 265 255 245 235 225 12:00 12:30 13:00 13:30 14:00 14:30 15:00 15:30 
Reflection Layer Height - km UTC on January 29, 2020 


Geometric Times Of Flight plotted in descending order 


and cropped to approximately match measured data. 14 


Measured 1-Hop Arrival Time is Consistent with 
Single Hop PHaRLAP Ray Trace Modeling 


PHaRLAP Simulation Data Provided by Nathaniel Frissell W2NAF 


TX: WWV 2020 Jan 29 12:00 UT RX: 1356.7 ke 
RX: WASFRF TX Origin: 40.68°N, -105.04°E; Azimuth: 153.8", Frequency: 5.000 MHz mye ese. 


Refraction height inferred from 10.0 
idealized geometry consistent 
PHaRLAP prediction. 


- GPSDO Timing + 


IC-7610 Proc Delay | 
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TX Power: 1.0 W, TX Gain: 1.0 dB, RX Gain: 1.0 dB 


Measured TOF 
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Ee cae Ts ae Predicted group delay near 5 mS consistent with measured arrival time of 4.93 mS 
Measured Data for primary 1-hop mode at 1201z. T = 9.68 (GPSDO) — 4.75 (IC-7610) = 4.93 mS. 
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Measured 2-Hop Arrival Time is Consistent with 
PHaRLAP Ray Trace Modeling 


PHaRLAP Simulation Data Provided by Nathaniel Frissell W2NAF 


TX: WWV 2020 Jan 29 13:24 UT 
RX: WA5FRF TX Origin: 40.68°N, -105.04°E; Azimuth: 153.8°, Frequency: 5.000 MHz 


RX: 1356.7 km 
Ray: 1355.7 km 
dx: 1.6 km 


ceiver 


Refraction height inferred from 
timing experiment consistent 
PHaRLAP prediction. 


Plasma Frequency [MHz] 
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960 .00000us 1.041 


Pi Le Piste [ecioe JP | Predicted group delay near 6 mS consistent with measured 2-hop arrival time of 
Measured Data 5.67 mS at 1329z. T = 9.46 (GPSDO) — 4.75 (IC-7610) + 0.96 (Extension) = 5.67 mS 
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Measured 3-Hop Arrival Time is Consistent with 


3F Proplab Pro 3D Ray Trace Modeling 


“E2- Agilent Technologies WED JAN 29 14:44:01 2020 
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The best measurement technique is to 
reference the leading cycle of the primary 
and last cycle of the delayed pulses since 
they are not distorted by superposition. 


RP 


WWV (40.0N/105.3W) to WASFRF (29.5N/98.9W), 1305 km, 5 MHz 
January 29, 2020 at 1430 UTC, from Proplab Pro V3 
ormnary wave, 3F mode, true pore circle Lisiaidl 151.5° 


Pe ee See tae ay 


: praca Pro V3 (IRI 2007) simulations 
- provided by Carl Luetzelschwab KOLA 
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Separation of WWV and WWVH 
During Simultaneous Reception 


¢ There are many times when the band is open to both WWYV in Colorado and WWYVH in Hawaii. 
Because they are on the same frequency the spectral contributions overlap in a spectrogram. 

¢ WWV and WWVH alternately transmit standard side tones of 500 and 600 Hz on most (but not all) 
minutes of each hour. See https://tf.nist.gov/stations/iform.html . Their use as a separation tool was 
recognized by David Kazdan AD8yY. 

e These signals are as precise as the carrier and can be used to separately measure Doppler data if 
recovered as an Upper Side Band (USB) tone. 


Minute WWV WWVH 
Even 500 Hz 600 Hz 
Odd 600 Hz 500 Hz 


Numerous separation schemes are possible using one or both tones. The option reported here is a fully 
deinterlaced technique that places 100% continuous data records for each station on separate outputs. 
The method requires two receivers with the ability to change frequency every minute so that one receiver 
continuously monitors WWV and the other WWVH. This method is also useful using only one 
programmable receiver to yield 100% data on just one station. 


MAIN = f. — 0.5; SUB = fe - 0.4; CHANGE every minute 


Deinterlace Scheme Using 


ANT 1 BW3.0k ) SFE 0) _/7\ 12:2036 ANT1 BW3.0k) SFT 0) (9 \ | 
. 1 4 3 5 7 BKIN $1 3 5 Lali 
Dual (or Two) Receivers to | ti 
4.999.50 
Separate WWV and WWVH re) saan) 


Even minute: 
MAIN on 500Hz for WWV 
SUB on 600Hz for WWVH 


Radio with A 1kHz Output Tone 


Dual Receivers A: always WWV 
(e.g., 1C-7610) 8B B: always WWVH 


4.999.50 translates pla 9 


Change A/B 500Hz USB to 1000Hz (EuaEsimmcia 1 RET USI | MARE | ERPOVEET 


Radio Interface 
(e.g., Cl-V) 
- Or Manual - 


4,999.60 translates ANT | MO eA 1221 a2 ANT! GK) Ge 
600Hz USB to 1000Hz 


BKIN 


USB FiL1 


4.999.60 


Rew) veBwW SPECTRUM SCOPE CENTER MAIN Grid 0.5k/10dB 
. . -2.5k +2.5k} -2.5k +2.5k 
Odd minute: 


DUAL.W’ 2) USB) FiL1 


Odd /Even e.g., LSB of Minute Digit 


Odd/Even Minute MAIN on 600Hz for WWV 
Determination 


Method can be 

(e.g., Radio automated if odd/even 
Controlled Clock) | Minutes can be sensed 
electronically. 


SUB on 500Hz for WWVH 


0 +1.0 +2:0 


N.B. Scheme also works with a single receiver with programmable VFO Eouleneay Camere iit al ili 
or memory channels to output either WWV or WWVH alone. IC-7610 screenshots in the two modes 


required to deinterlace the WWV/WWYVH side 
tones. The radio was stabilized with a GPSDO. 


Completely Deinterlaced WWV and WWVH Data 
- |IC-7610 Dual Receivers Using CHANGE to Manually Switch MAIN and SUB — 
- Data is Present for Each Tone Every Minute It is On the Air - 


Only WWV on MAIN Only WWVH on SUB 
(Multiple Modes Present) (Single Mode Present) 


Conclusions and Acknowledgements 


Using experimental tools targeting specific ionospheric behavior can help 
decipher the multiple physical processes involved. The purpose of this paper was 
to give examples of some of the tools developed for the HamSCl community. 


The results of the studies conducted to date support that at least some of the 
observed mode splitting may result from different Doppler shifts for multiple hop 
modes. But the data shows much more complexity than can be accounted for by 
this one mechanism. Alternate explanations have been proposed and are the 
subject of future research. An automated TOF data collection system with 
simultaneous spectrogram data would be of immense help. 
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